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Estimating obstruction fading on line-of-sight microwave radio 
paths requires statistics on positive refractivity gradients in the 
lowest 100 meters of the atmosphere. In this paper, we describe a 
semiempirical climatological model that predicts quantitatively the 
occurrence frequency of positive refractivity gradients and the geo- 
graphical variation of this frequency within the contiguous United 
States. The occurrence frequency is parameterized by six factors, for 
which numerical values are obtained from available climatological 
and physiographical data. By modeling the two diurnal atmospheric 
boundary layer regimes with separate probability density functions, 
we produce cumulative probability distributions for refractivity gra- 
dients. The model is normalized using measured refractivity gradient 
distributions available for 17 sites within the United States. The 
model can be used for line-of-sight microwave radio path engineer- 
ing; a companion paper presents a method for determining antenna 
heights. 

I. INTRODUCTION 

Line-of-sight microwave radio transmission performance is suscep- 
tible to atmospheric conditions, in particular to spatial variations in 
the microwave refractive index. Nearly always during the day and 
frequently at night, sufficient vertical mixing occurs of the low-level 
air layers, through which the microwave beam propagates, so that 
vertical refractive index gradients encountered are small. However, 
horizontal layers containing strong vertical gradients in air tempera- 
ture and water vapor pressure, which determine the refractive index, 
sometimes develop at night. The presence of these layers can cause 
the received signal to fade (decrease in strength). 

This paper addresses the meteorology of obstruction fading, which 
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is caused by a large positive gradient in the refractive index, which is, 
in turn, caused by a large positive water vapor pressure gradient. The 
ray path between the transmitting and the receiving antennas is lower 
than normal in the presence of a positive refractive index gradient, 
and the terrain can block the ray path when the positive gradient is 
sufficiently large. Although infrequent, the consequences of obstruc- 
tion fading are severe: transmitted signal is steadily blocked for as long 
as an hour at a time. Obstruction fading can be avoided by locating 
microwave antennas sufficiently high above the ground. 

Prolonged fades can also be caused by strong negative gradients that 
refract the transmitted microwave beam downward to such an extent 
that it ceases to illuminate the receiving antenna. Climatic conditions 
and path geometries in the United States are such that this is usually 
a secondary phenomenon compared to obstruction fading. 

The presence of negative gradients of more modest proportions is 
associated with multiple propagation paths through the atmosphere, 
which cause phase interference fading (multipath fading). Phase inter- 
ference fading is a major factor in microwave system design. The 
meteorology of multipath fading is a separate topic, and is not ad- 
dressed in this paper. 

Probability distributions for the occurrence of refractive index gra- 
dients have been compiled by Samson, 1 Bean et al., 2 and Dougherty 3 
from measurements of vertical profiles of meteorological parameters. 
The distributions are intended for use in engineering microwave paths. 
However, these sets of compilations are for isolated locations, and they 
have not been generalized to account for the varied climates to which 
different microwave routes are exposed. 

The purpose of the present work is to devise a climatological model 
to predict quantitatively the geographical variation of the occurrence 
frequency of positive refractive index gradients within the contiguous 
United States. The results can be used to estimate obstruction fading 
occurrence on microwave radio paths to help determine optimum 
microwave antenna heights. The result of the model is a predicted 
cumulative probability distribution for positive refractive index gra- 
dients calculated for an arbitrary location within the contiguous United 
States for arbitrarily large gradients. The present model represents the 
first attempt to predict refractivity gradient probability distributions 
from climatological data. Ikegami et al. 4 and Akiyama and Sasaki 5 
have performed fits to measured refractive index gradient distributions 
but without meteorological parameterization. 

This model is designed to provide predictions for the 48 contiguous 
United States. Two submodels must be used to account for the varied 
climates experienced by this large geographical entity: a coastal region 
submodel and an interior region submodel. The interior region sub- 
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model is based on a relatively small experimental data base and the 
reliability of the overall model for interior regions is possibly lower 
than for coastal regions. However, since obstruction fading occurs 
mostly in coastal regions, the model is most reliable where it is most 
useful. 

The predictions of this model are compared to measured distribu- 
tions 1 and, in general, show good agreement. We attribute the imper- 
fect agreement that exists in a few cases to local effects that cannot be 
accounted for without introducing additional complexity with the 
attendant risk of making the model unwieldly to implement. That is, 
we use a geographical resolution for the model, which has been judged 
optimum by weighing resolution improvements against simplicity of 
model use. 

II. MODEL FORMULATION 

The model is formulated in two stages. First, we choose parameters 
(Section 2.1) based on meteorological principles, to relate the occur- 
rence frequency and magnitude of refractive index gradients to avail- 
able climatological and physiographical data. Second, we develop a 
method for producing a refractive index gradient probability distribu- 
tion from these parameters (Section 2.2). The model is quantified by 
comparison with measured distributions. 1 We make final minor ad- 
justments to the model to provide consistency with information on the 
geographical distribution of obstruction fading gleaned from transmis- 
sion performance experience in the Bell System. 6 

2.1 Meteorology 

The refractive index, n, for microwaves can be expressed as refrac- 
tivity, 7 

N = 10 6 (n - 1) = ^ (p + 4810 j\ , (1) 

where N is the refractivity in N-units, T is the air temperature in 
kelvin, P is the air pressure in millibars, and e is the water vapor 
pressure in millibars. Water vapor pressure is the governing contributor 
to producing positive refractivity gradients (A/') at night. Therefore, a 
positive refractivity gradient model must predict climatologically the 
frequency of occurrence of positive e gradients. 

During obstruction fading, the important part of the microwave 
beam generally lies within the lowest 100 m of the atmosphere and 
available radio meteorological data have been obtained, therefore, for 
100 m. However, since probability distributions are functions of height 
increments, we present a method (Section 2.2) for using a height 
increment appropriate to the particular application. 
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Several mechanisms can produce inverted (positive gradient) water 
vapor vertical profiles: 

(i) Condensation of water vapor at the bottom of a surface-based 
nocturnal temperature inversion, e.g., dew or ground fog, causes ab- 
solute humidity to decrease there. 8 However, the amount of water 
vapor condensed is insufficient to produce large gradients across a 
layer as thick as 100 m. 

(ii) Strong convection processes in relatively dry air can produce 
water vapor inversions within the daytime boundary layer. 9 Develop- 
ment of a nocturnal temperature inversion could cause the water vapor 
inversion to persist throughout some of the night. The magnitude of 
these gradients is insufficient to cause obstruction fading, however. 

(Hi) Uplands generally receive more rainfall than surrounding low- 
lands, provided the altitude difference is large enough. When situated 
within relatively dry lowland regions, daytime evaporation from the 
uplands can produce elevated moist air layers. Such an effect has been 
observed. 10 Nocturnal radiative temperature inversions could cause 
the water vapor inversion to persist during the night. Such a mecha- 
nism is greatly limited geographically, however. 

(iv) Nocturnal advection of moist air can produce inverted water 
vapor profiles. Nocturnal land breezes at water body-land interfaces, 
are a possible mechanism. However, the effects of land breezes usually 
do not penetrate more than a few tens of kilometers inland and, thus, 
cannot explain severe obstruction fading observed farther inland. 

(v) We propose large-scale advection of moist air as the mechanism 
for producing most of the obstruction fading on line-of-sight microwave 
radio routes. This mechanism is supported by various characteristics 
of the limited information on observed obstruction fading in the Bell 
System: 6 (1) Often large areas are affected by a single episode. (2) The 
length of microwave paths requires vertical gradients to be horizontally 
homogeneous over approximately 40-km distances to produce fading. 
(3) Relatively flat areas are especially susceptible; flat areas are con- 
ducive to large-scale, low-level advection at night. (4) Regions within 
300 km of coastlines are especially susceptible; 300 km is the distance 
air could travel at typical speeds at approximately 100-m altitude 
during 12 night hours. (5) Finally, a study of a severe fading episode in 
the Bell System showed that the onset of obstruction fading occurred 
later for greater distances from the coast. 

An understanding of the details of the proposed large-scale moisture 
advection mechanism requires knowledge of atmospheric motion on 
various distance scales. We will now discuss atmospheric motion within 
the boundary layer, across the continent, and within air masses. 

Obstruction fading usually occurs during fine weather conditions, 
when anticyclonic air masses envelop the affected region. Thus, the 

806 THE BELL SYSTEM TECHNICAL JOURNAL, JULY-AUGUST 1 981 



characteristics of the boundary layer and their diurnal variation only 
in fine weather conditions need be considered. 

During the day, solar radiation causes convective mixing which 
produces a relatively deep (==1 km) boundary layer which is vertically 
well mixed throughout. Thus, vertical gradients of water vapor and 
horizontal momentum are relatively small across the lowest 200 m 
where the microwave beam lies. Small vertical gradients of horizontal 
momentum produces small wind speed and direction gradients (Fig. 

la). 

At night, radiative cooling of the earth's surface produces a relatively 
thin (as thin as ~100 m) density-stable boundary layer. The high 
stability suppresses turbulence and very little vertical mixing occurs; 
water vapor inversions would be maintained once produced. Relatively 
large wind speed and direction gradients could be maintained across 
the boundary layer 11 (Fig. lb) because of the limited vertical mixing of 
horizontal momentum. Calm winds often occur at the surface. 

Large-scale air motion is determined by the midtropospheric (~6 
km) winds. These currents steer air masses across the North American 
continent, in general, from west to east. Air masses that affect the 
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Fig. 1 — Air flow in lowest 200 m. Arrows indicate relative wind speed and directions 
at two heights, illustrating the enhancement of differential flow at night which contrib- 
utes to moisture inversions in coastal regions. 
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United States usually originate over northwest Canada or the eastern 
Pacific Ocean. 12 Pacific air masses crossing the continental divide lose 
much of their moisture through lifting-induced precipitation over the 
western mountains. Therefore, except for the west coast, air masses 
are relatively dry until they reach the Great Lakes, Atlantic Ocean, or 
Gulf of Mexico. Once they reach these large water bodies, evaporation 
increases their moisture content. 

Air flow at moderately low levels (:S2 km) within anticyclonic air 
masses is essentially circular (for the purpose of this discussion) with 
respect to the air mass core, and in a clockwise direction. Thus, 
moderately low level winds at a given location change direction as an 
air mass traverses the region. Furthermore, winds near the center of 
the air mass are very light because of small horizontal pressure 
gradients; therefore, highly stratified air layers can readily form there. 

The mechanism for producing water vapor inversions by the large- 
scale moisture advection that we are proposing is as follows: dry air 
masses approach the east (south) coast of the United States from the 
west (northwest). Before the air masses reach the coasts, moderately 
low level (:S2 km) wind flow is mainly offshore in coastal regions, thus, 
maintaining dry air there. As the air mass core comes sufficiently close 
to the coast, winds to the south and west of the air mass core shift to 
an onshore direction, advecting moist air inland. If the offshore- 
onshore wind shift occurs near sunset, the characteristics of the night- 
time boundary layer discussed earlier, wind speed and direction 
differences across it, cause moist air to be advected over the top of 
surface-based dry air. The differential wind speed advects moist air 
inland faster at approximately 200 m than at the surface. The dif- 
ferential wind direction, which is normally a clockwise shift with 
height, also contributes to maintaining faster inland penetration of 
moist air at approximately 200 m relative to the surface. 

Production of positive e gradients at a particular site requires the 
site to be sufficiently close to a moisture source so that moist air can 
reach the site during the nocturnal stratified air regime. Measured 
refractivity gradient distributions 1-3 show that elevated moist layers 
occur most frequently near large bodies of water such as oceans, gulfs, 
and large lakes. However, measurements show also that elevated moist 
layers occur at sites too far from these large water bodies to have 
had them for their source of moisture. Thus, mechanisms other than 
large-scale advection probably also contribute to positive refractivity 
gradients. Each mechanism cannot be modeled specifically, as for 
large-scale moisture advection. Therefore, a large-scale moisture-ad- 
vection model is used near large water bodies (coastal regions), while 
a mechanism-independent model is used for other areas (interior 
regions). 
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2.1.1 Coastal region 

The coastal region model is parameterized by six factors relating the 
occurrence of positive e gradients to parameters obtainable from 
seasonal climatological and physiographical maps. The factors are 
frequency of temperature inversion occurrence, fr, air homogeneity, H; 
water body proximity, D; frequency of wind direction occurrence, f w ; 
air moisture capacity, V; and surface moisture availability, M. The 
proportionality among the six factors and their individual normaliza- 
tion is obtained empirically from measured N' probability distributions 
(Section 2.2). J The factors /i and H specify the occurrence of stratified 
air, D and f w specify the production of positive e gradients, and V and 
M specify the magnitude of the positive e gradients. 

The occurrence of stratified air in the lowest 100 m coincides with 
the existence of temperature inversions (temperature increase with 
height). Hosier has compiled maps of the occurrence frequency of 
surface-based temperature inversions for four seasons within the 
United States. 13 The frequency of temperature inversion occurrence, 
fi, whose units represent the fraction of total time that inversions occur 
for a given location and season, is employed as a factor in the coastal 
model. 

For obstruction fading to occur, large positive e gradients must occur 
over the entire microwave path (a standard path length of 40 km is 
assumed in this work). Thus, the horizontal homogeneity of the strat- 
ified air, which is determined by the flatness of the terrain profile, 
must be large for the occurrence of obstruction fading. Therefore, a 
factor, H, is introduced that represents the horizontal homogeneity of 
the stratified air. Values for this factor are obtained from a map in the 
United States National Atlas, "Classes of Land Surfaces Forms." 14 In 
this map, each location is represented by three characteristics: slope, 
local relief, and profile type. In the present model, each characteristic 
is defined as a subfactor of H, 

H = HxHJIz. (2) 

Numerical values of Hi, H 2 , and H 3 are acquired, after all other 
empirically determined parameters are obtained, by comparison of the 
model with the geographical distribution of microwave transmission 
performance in the Bell System. 6 Transmission performance experi- 
ence is employed since it represents the only existing path- averaged 
obstruction fading information. Table I gives the values obtained for 
H u H 2 , and H 3 . 

It is useful to compare refractivity gradient distributions which are 
predicted by the present model and represent path averages, with 
measured distributions compiled at single points. In such cases path 
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Table I — Air homogeneity subfactors. 
Classes are obtained from Ref . 1 4. 
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averaging is removed from the predictions. This is accomplished by 
setting H to its maximum value, 1.0, in those instances. 

For the coastal region, production of positive e gradients requires 
proximity to a large water body. Moist air originating at the coastline 
at sunset can penetrate a maximum distance inland during the noc- 
turnal stratified air regime. Therefore, we propose the following 
smoothed-step distance function to represent this effect: 

D = [1 + (d/do) 4 T\ (3) 

where d is the minimum distance to an ocean, gulf, or one of the Great 
Lakes and d = 300 km, as determined by the maximum distance 
inland that the air can penetrate during 12 hours when moving at an 
average speed at a height of approximately 100 m. 

For every location there is a wind direction that maximizes the 
relative advection of moist air over dry air. The appropriate direction 
varies with distance to and shape of the coastline nearest the site. The 
procedure for specifying this direction for different sites is discussed in 
Appendix A. We propose that the production of positive e gradients is 
proportional to the frequency, f w , that wind blows along this direction. 
Values for f w , in units of fraction of total time, are obtained from 
frequency distributions of surface wind direction. 

The magnitude of the e gradients formed depends upon the water 
vapor capacity of the air. The water vapor capacity, in turn, depends 
upon air temperature, T. The saturation vapor pressure function, 
e s (T), is a logical choice to represent the dependence of water vapor 
capacity on air temperature. It is proposed that a water vapor capacity 
factor be represented by 

V = , (4a) 

e 8 o 

where e s (T) is in millibars and e s0 = 50 mb, an approximate upper 
limit for naturally occurring vapor pressures. The function e s (T) can 
be represented by 8 

e s (T) = exp[21.71 - 5433/T], (4b) 
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where T is in kelvin. Temperatures used are those obtained from 
monthly average temperature maps of the United States. 

The magnitude of the e gradients produced by advection from large 
water bodies is reduced if substantial smaller local sources of water 
exist at the site. This is caused by surface water evaporating into the 
dry lowest air layer. This effect is modeled with a surface moisture 
factor 

Mc ^-, (5) 

mo + m 

where m is a surface moisture index obtained from a map of moisture 
regions, 17 m = 70, and mi = 20, obtained by comparison with measured 
N' distributions. 1 The subscript C denotes the coastal region since a 
different moisture factor is employed for the interior region. 

In summary, the six factors used to determine the occurrence of 
positive e gradients in the coastal region are ft, H, D, f w , V, and M c . 

2.1.2 Interior region 

The model for interior sites is a modification of that developed for 
the coastal sites. The two factors representing the production of 
stratified air layers, H and ft, are the same as for the coastal model. 
However, the two factors representing the production of positive e 
gradients, D and f w , are combined into a single constant factor, A'. 
This is required because it is not feasible to take into account the 
variety of mechanisms thought to contribute to gradient production in 
interior regions. Determination of A' is discussed below. Finally, of the 
two factors representing the magnitude of the e gradient produced, V 
and M, only M need be modified. This modification is required since 
the surface moisture index, m, provides the necessary information 
regarding the availability of moisture from local sources. Thus, the 
proportionality between M and m is inverted with respect to the 
coastal site model. The interior region model for M is 

Mi = . (6a) 

rn.3 

Since the relative advection factor, A', is constant, m^ can be absorbed 
into it yielding a = A'/m 3 . We determine empirically that a = 4.7 X 
10 -5 and m 2 = 100 (Appendix B). Equation (6a) can be rewritten as 

Mi = m + m.2. (6b) 

The five factors used to determine the occurrence of positive e 
gradients in the interior region are ft, H, a, V, and Mi. 

The interior model, instead of the coastal model, is applied to the 
Pacific Coast of the United States. The large-scale moisture-advection 
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mechanism upon which the coastal model is based is not valid for the 
Pacific Coast since air masses which cross the coast already contain 
much moisture. Furthermore, relatively high coastal mountain ranges 
prevent low-level nocturnal differential advection. 

2.2 Re fr activity gradient probability distributions 

Examination of measured N' cumulative probability distributions 1 " 3 
indicates that the shape of each curve is characterized by three features 
when plotted on normal probability axes (Fig. 2): a central linear 
segment and positive and negative tails. Ikegami et al. 4 and Akiyama 
and Sasaki 5 have noted these features and have attempted to fit each 
with different mathematical formulas, but no attempt was made to 
interpret them meteorologically. We hypothesize that the central 
segment is produced by an atmospheric regime characterized by at 
least some vertical air mixing, and the tails are produced by a regime 
characterized by stratified air. In general, the tails are distributed 
asymmetrically about the mean. This is expected, since different 
meteorological mechanisms produce positive and negative gradients in 
stratified air. Therefore, a mixture of at least three distributions is 
required to represent the total cumulative probability distribution 
completely. In view of the complexity of handling a mixture of three 
distributions relative to two distributions, and since only positive N' 
gradients are of concern for obstruction fading, we do not attempt to 
model that half of the cumulative probability distribution where N' is 
more negative than the mean. Therefore, we employ a mixture of two 
distributions. 
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Fig. 2 — One of 17 measured refractivity gradient distributions compiled by Samson. 1 
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In general, the mixed regime distribution comprises approximately 
80 percent of the total distribution (see discussion of p* below). 
Furthermore, it is well documented that mean refractivity gradients in 
mixed air are similar to those expected from normal daytime vertical 
distributions of water vapor, temperature, and air pressure. 718 There- 
fore, since gradients developed during the stratified regime evolve from 
those existing during the mixed regime, and since the negative-going 
half of the total distribution is not modeled (relieving the need to 
account for asymmetry), it is reasonable to assume that each of the 
two component distributions can be represented by normal distribu- 
tions that share the same mean, ju. Figure 3 illustrates this assumption. 
These component distributions are given by the probability density 
functions 



9 xn r ifN'-u 

p m (N') - (2tt <x^)- 1/2 exp - - £ 

. 2_ 

1 /AT' _ „\ 

p 8 (N') = (2tt a 2 r 1/2 exp 



l/N'-u 



(N'>ii), (7a) 



(AT=>/i>. (7b) 



2 \ a s 

where m and s denote the mixed and stratified regimes, respectively. 
The total probability density function is given by 

p(N') = (l-p*)Pm+p*p s , (8) 

where p* is the proportion of mixture for the two component distri- 
butions. 

We empirically compared a mixture of two distributions with Sam- 
son's U.S. data and determined thatp* = 0.2 and a m = 15 TV-units/km 
and that both are essentially invariant with respect to location and 
season. Bean et al. have compiled seasonal maps of the average 
refractivity gradient in the lowest 1-km depth. 2 This data is used to 
specify /i as a function of location. 

The standard deviation for the stratified regime, o s , is parameterized 
separately for the coastal and interior regions by 

o sC = a (fiH) i/2 (Df w VM c y /4 (9a) 

and 

ff./=ao(/iff) 1/2 (aVM/) l/4 , (9b) 

where ao is a constant for a given layer thickness and path length, and 
the exponents, M> and V*, are determined empirically. Both o»c and o 8 i 
are calculated for each site and season, and the larger of the two is 
employed to calculate p s for each season. Total probability, p(N'), is 
calculated for each season and the seasonal probabilities are averaged 
to produce an annual probability density function for each site. Finally, 
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Fig. 3 — Illustration of two probability density functions employed to represent a 
refractivity gradient distribution. Dashed portions of the functions denote halves which 
are not modeled. The dotted line denotes a possible "true" negative-going half of the 
stratified distribution which is, in general, asymmetric about the mean. 

the probability density functions are integrated from +oo to //. to 
produce N' cumulative probability distributions for positive N'. 

Samson's 17 measured N' distributions 1 were used to determine 
empirically that a for a 100-m-deep layer is 440 iV-units/km. N' 
distributions have been compiled at several sites for surface layer 
thicknesses (Az) different from 100 m. Samson compiled annual dis- 
tributions at Cardington, England for 75-, 150-, and 500-m- thick layers. 1 
Dougherty compiled annual distributions at Cocoa Beach, Florida for 
50- and 100-m-thick layers. 3 Values for the climate parameters appro- 
priate to each site were used in an empirical determination of oo for 
the five available layer thicknesses. H=l was used; the reason for this 
choice has been discussed in Section 2.1.1. 

The results are plotted against (Az)~ 1/4 in Fig. 4 and indicate such a 
relationship. The line in Fig. 4 represents cto(Az) = 1400/Az 1/4 , where 
Az is in meters and ao is in iV-units/km. Thus, o*o for an arbitrary layer 
thickness, can be obtained from this plot. 

III. RESULTS AND DISCUSSION 

Complete specification of a predicted N' probability distribution is 
made by the parameters p*, ju, a m , and ov The parameters p* and a m 
are constants; however, ju and a s must be specified for each season as 
a function of location. Bean et al. provide seasonal maps for ju mea- 
surements, 2 and Table II summarizes specification of ov 

It is useful to compare a few predicted distributions with measured 
ones to check the performance of the model. Figures 5 through 7 show 
examples of N' cumulative probability distributions predicted by the 
model for 100-m layers (thicknesses for which measured distributions 
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Fig. 4— Dependence of a upon surface layer thickness, A 2. Data points are derived 
from Samson s compilations' for Cocoa Beach, Florida (•), and Cardington, England 
(O). 

are available); measured curves are also shown. Predicted curves were 
calculated using H = 1 (see Section 2.1.1). The ability of the model to 
successfully predict large differences in distribution shapes for two 
sites separated by only 500 km is illustrated by comparing the Cocoa 
Beach and Key West curves. The model's ability to respond to seasonal 
variations in both coastal and interior regions is illustrated for Cocoa 
Beach and Denver. 

The agreement between predictions and measurements is generally 
good; however, it is not perfect. Each measured distribution, with a 
few exceptions, was compiled from approximately 900 measurements 
(Table III). Thus, the data point density on the distribution tails is 

Table II — Summary of method for evaluating a s parameters. 
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Fig. 5 — Comparison of positive refractivity gradient distributions predicted in this 
work (solid line) with those compiled by Samson 1 from measurements (dashed line) for 
a 100-m-deep layer for coastal site Cocoa Beach. 

low. Therefore, some contribution to the disagreement is expected 
from statistical noise. 

A more complete assessment of the agreement between predictions 
and measurements has been made by plotting predicted vs measured 
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Fig. 6 — Comparison of positive refractivity gradient distributions predicted in this 
work (solid line) with those compiled by Samson 1 from measurements (dashed line) for 
a 100-m-deep layer for coastal site Key West. 
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Fig. 7— Comparison of positive refractivity gradient distributions predicted in this 
work (solid line) with those compiled by Samson' from measurements (dashed line) for 
a 100-m-deep layer for Denver (an interior site). 

annual probability for N' 2: 100 N-units/km for Samson's 17 sites (Fig. 
8). 1 Such an evaluation requires that path averaging be removed from 
the predictions for comparison with point measurements. This is 
accomplished by setting the air horizontal homogeneity factor, H, to 
its maximum value (1.0) for all sites. Therefore, probabilities calculated 
in this manner are characteristic of a single point along the 40-km path 

length. 

It is most useful to compare predicted and measured distributions 
for the largest refractivity gradients for which measured probabilities 
could be obtained from Samson's data. 1 As the optimum gradient for 
comparison, 100 N-units/km was selected. However, not all of the 
measured distributions for each season at a given site intersected 100 
N-units/km. Therefore, for those that did not, distribution tails were 
extrapolated linearly on the original graphs to estimate probabilities 
at 100 N-units/km for those seasons. Annual values of P(N' > 100 N- 
units/km) were obtained by averaging the four seasonal probabilities. 

Poorer agreement between predicted and measured distributions is 
expected for interior sites, since a single specific mechanism has not 
been modeled as for the coastal sites. However, obstruction fading is 
not a serious problem in interior regions. For all sites except four, the 
prediction is near the measurement. Two of the four sites are interior 
ones, discrepancies for which were anticipated. The four anomalous 
sites (Denver, El Paso, New York, and Long Beach) all show consid- 
erably higher occurrence frequency of positive gradients than pre- 
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Table III — Sites for which Samson's measured refractivity gradient 

cumulative probability distributions 1 were employed and approximate 

number of individual observations associated with each. Symbols in 

last column identify points in Figure 6. 
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Brownsville, TX 
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Miami, FL 
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Charleston, SC 
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New York, NY 
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Cocoa Beach, FL 
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880 


N 


Denver, CO 


880 


F 
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Juliet, IL 
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Washington, DC 


840 
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dieted. This observation suggests that local meteorological effects not 
incorporated in the model are contributing to producing gradients. 
Indeed, all four measurement sites are unusual. The Denver site is on 
relatively flat ground near the foothills of the Rocky Mountains; the 
El Paso site is on a mesa approximately 60 m above flat, irrigated land; 
the New York site is located in a complex metropolis-ocean-land-bay 
region; and the Long Beach site is subjected to Santa Ana winds. 19 
Therefore, the anomaly is attributed to local effects on the probability 
distributions measured at a single point. If it were possible to measure 
path-averaged distributions for the four anomalous sites, the model 
predictions for them might, indeed, show better agreement. 

IV. SUMMARY 

This work presents a climatological model for arbitrary locations in 
the United States. It provides refractivity gradient occurrence distri- 
butions for positive gradients. The results fill in the gaps among the 17 
locations for which measured distributions have been compiled. The 
significance of this additional detail is that climate variations can be 
accounted for with much better geographic resolution than had pre- 
viously been possible, in the specification of microwave antenna 
heights. A method which uses these results to determine antenna 
heights has been developed by Vigants and is described in a companion 
paper. 20 
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Fig. 8 — Comparison of model predictions with Samson's measurements for 17 sites. 1 
Values plotted are cumulative probabilities for N' > 100 AT-units/km. The line represents 
perfect agreement between predictions and measurements. Sites are denoted by symbols 
which are identified in Table III. The four anomalously low points are attributed to the 
inability of the model to account for unusual local effects. 



APPENDIX A 

Critical Wind Direction Choice 

Climatological information does not exist on the occurrence fre- 
quency of wind shifts from continental to maritime directions as air 
masses traverse the coastlines. However, there is a certain wind 
direction relative to a coastline which should be related to the occur- 
rence of large-scale moisture advection. 

Before air masses reach the coasts, moderately low-level air flow is 
roughly parallel to the coast with land to the right of the wind vector. 
Furthermore, this direction is conducive to producing moisture gra- 
dients by the mechanism of wind direction variation with height, as 
discussed in Section 2.1. Therefore, the surface wind direction that 
most likely causes flow of marine air over continental air adjacent to 
the coastline is the direction parallel to the coast with land to the right 
of the wind vector. This is the critical wind direction at the coastline 
for an ideal linear coastline of infinite extent. In reality, however, the 
critical direction choice must account for the curvature of the coastline. 
This is required to ensure that air assumed to be continental (mari- 
time) in origin indeed is underlain by land (water) for a significant 
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distance upwind from the site. An appropriate choice of minimum 
distance is the distance travelled by air in one day at an average speed, 
approximately 300 km. 

Because of unavoidable subjectivity required in selecting critical 
wind directions for some sites, the selection technique is best illustrated 
with examples. Figure 9 shows nine sites along the southeast U.S. 
coast, three each along lines drawn perpendicular to the Atlantic 
Ocean coast from three different coastal sites: A, B, and C. Critical 
directions for line AA' are parallel to the coast because land lies 
upwind along this vector for a distance of at least 300 km. This assures 
that air arriving at sites along line AA' is indeed continental. Such an 
assurance is not obtained for sites along line BB' if the critical 
directions are chosen parallel to the coast at B. The directions are 
chosen so that there is a 300-km air flow trajectory flow over land 
upwind from sites along BB', as shown in Fig. 9. 

The situation at site C is complementary to that at site B. Here, the 
critical wind direction is chosen so that the wind vector is underlain 
by water for a distance of 300 km upwind from site C. The critical 
direction at the other two sites along CC is chosen assuming that 
maritime air has its source over the Atlantic Ocean rather than the 
Gulf of Mexico. This is done since a gulf source would require a critical 
direction at site C, e.g., as shown by the dotted arrow in Fig. 9. This 
direction, in turn, would imply continental air originating in the 
southern half of the Florida peninsula. Obviously, this cannot be 
considered a source of dry, continental air. Critical direction choices 
illustrated in the above examples are typical of those made for nearly 
all locations of the United States. 

APPENDIX B 

Empirical Determination of a and m 2 

The values of a and m.2 for the interior regions are obtained by 
comparing measured refractivity gradient probability distributions for 
Columbia, Missouri. 1 The only other interior sites for which measured 
distributions are available (Denver and El Paso) are characterized by 
unusual physiographical features (Section III) and are deemed unsuit- 
able for specifying a and m.2. Values are obtained for fi, ft, V, and m for 
Columbia. H = 1 is used for comparing the model with point gradient 
measurements (see Section 2.1.1). The value for ao obtained for coastal 
regions is used. From eq. (9b), a s / for Columbia is given by o 8 i = 
constant X a (30 + 7712), where the constant is evaluated from the above 
values. 

Determining the optimum agreement between calculated N' distri- 
butions (using ct s / as specified above) and the measured N' distribution 
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Fig. 9— Examples of critical wind-direction choice. Dashed lines are perpendicular to 
e Atlantic Coast. Arrows represent critical wind directions chosen for sites at their 



the Atlantic Coast. Arrows represent 

tips. Arrow length is 300 km (see Appendix A). 



for Columbia, yields a relationship between a and m 2 . We determine 
the remaining variable by subjecting calculated values of a s i{m) to 
constraints for the possible range of m: -50 < m < 120. Columbia's 
surface moisture index 17 of 30 is intermediate within this range. A 
second point to completely specify a and m 2 is obtained by assuming 
that a s i(m = 120) for Columbia is equal to the smallest value of a sC (jn 
= 120) calculated for all of the 14 noninterior sites. Thus, we determine 
that a - 4.7 X 10" 5 and ra 2 = 100. 
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